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Abstract 



A simple macroscopic continuum elasticity theory (CET) is used to calculate 
the structure relaxation of single- wall carbon nanotube (SWNT), an analytic 
formula is obtained. We also expand an atomic scale three-parameter empir- 
ical model [ T. Lenosky et al. Nature 355, 333(1992)] in order to correctly 
describe the bond-length change effects. The structure relaxation of SWNT 
expected by the model is good in agreement with our CET results, and very 
well consistent with the previous calculation from a first principles local den- 
sity function approximation. Using the expanded Lenosky model, we calculate 
the strain energy of bending tube. The obtained results are good in agree- 
ment with the previous theoretical expectation. It shows the model may be a 
good simple replacement of some more sophisticated methods on determining 
carbon networks deformations. 
PACS : 61.46. +w, 61.48. +c, 63.20.Dj 
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Recently, carbon nanotubes with cylindrical graphite structures have been intensively in- 
vestigated JI|, and many experimental and theoretical researches have been performed 
The structure of nanotubes is qualitatively well known through the simple construction of 
rolling a perfect graphene sheet, where only one parameter is to be determined: The lat- 
tice parameter or a bond-length. The first-principles local density function approximation 
(LDA) |||7]], tight-binding (TB) approximation [|] and empirical potential methods [^-|TT| 
have been used to determine the bond-length of tube. In all these theoretical studies, a small 
relaxation effect has been found, the bond length of tubes is not equal to that of graphene 
sheet since the curvature of the tubes and thus the structural inequivalence between the 
axis and normal direction render the carbon-atom hexagons distorted, although the details 
of the bond relaxation from different authors are different. In this letter, we use a sim- 
ple macroscopic continuum elasticity theory (CET) and an expand microscopic empirical 
model to calculate the bond length relaxation of single- wall carbon nanotubes (SWNTs). 
We find from both the two methods: (1) the radius of any SWNT is slightly larger than 
the expected value from the rolling graphene sheet. (2) The radius relative increasing of 
SWNT is reverse proportion to the square of the tube radius, and independent or slightly 
dependent on the helicity of tubes. Our results are in good in agreement with the previous 
LDA results. (3) Although the energy contribution of bond- length change is very small in 
straight SWNT Ul2"f , it can strongly affect the structure relaxation of tubes. (4) We use the 
expand model to describe the bending of nanotube, and find the obtained results are good 
consistent with the previous theoretical expectation. 

SWNT can be indexed by a pair of integers (n 1; n 2 ) 0. From the rolling graphene 
model, we have 



a/ (n\ + rix ■ n 2 + nl) ■ a 
Po ~ 2^ ' (1) 

where p is the non-relaxed radius of SWNT, a = V3r is the lattice constant of graphite 
sheet, r = 1.42 A is the bond length of graphene. Supposing the relaxed radius is p, we 
define 

e = ^. (2) 
Po 

A number of theoretical and experimental studies show the deformations of graphenes 
and SWNTs can be well described by the CET ||7|,p|, p!0| , p!3| , ^4]] . The rolling energy of graphene 



is only dependent on the rolling radius, and the stretching energy is only dependent on the 
relative elongation, independent on the rolling and stretching direction of graphite sheet, 
and we have 

E r = C/p 2 , (3) 
Es = \E" ■ e\ (4) 



Here E r and E s are the rolling and stretching energy per atom, respectively. C ~ 1.2 ~ 
2.0 eVA 2 , depending on the different models, E" « 58 ~ 59 eV 0,|,|l|Jl|,[ll . When rolling 
a graphite sheet to SWNT, the bonds are stretched in the circumference direction of tube, 
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(6) 



( we neglect the possible smaller strains in the axis direction). The strain energy per atom 
is the sums of the rolling and stretching energy. 

*-^ + ^ < 5 > 

Minimizing the strain energy, we obtain 

2C 1 2C 1 

So we obtain, e ~ (selecting C = 1.44 eVA 2 JTJJ )where the unit of p is ^- The results 
show: (1) any tube is relaxed when it is rolled from graphite sheet and the relaxation is slight 
but obvious; It is consistent with the previous theoretical expection fl,|5|,fryi~5| (2) Although 
the stretching energy from the dilation of tube radius E s is far smaller than the curvature 
energy contribution E r , can be safely neglected [12,lJ], the size of the radius dilation e is 
sensitively dependent on E s . 

For more carefully discuss the relaxation effects, we also consider a microscopic model. 
Lenosky et al. employed a three-parameter empirical model which can well describe the 
structure of the curved carbon networks while the bond lengths are constrained to r . Ou- 



Yang et al. |13| reduced the model to a two-elastic-modulus macroscopic continuous model 
and well treated the deformed energy of multi-wall carbon nanotube. However, for describ- 
ing some deformations which the bond lengths are obvious changed, such as bending or 
stretching tube, we need a more complete model which must including the contribution 



from the bond length [12], 



to') * U) 

{hi ■ Uij)(n,j • Uij) (7) 

(y) fe") 

where unit vector pointing from carbon atom % to its neighbor j, is the distance 

between the atom % and atom j, and hi is a unit vector normal to the fullerene surface at 
atom i. The summation ~Ylu\ is taken over the three nearest neighbor j atoms on atom i, 
and the summation ^fa) i s taken over only the nearest neighbor atom pairs. Supposing 
r ij — r o? fitting to LDA calculation, Lenosky obtained the values of (Si, S 2 , S3) which are 
(0.96,1.29,0.05) eV, respectively p|. 



Using an empirical 4th-neighbor force-constant model of graphite ||16|| , which has been 
well described the different deformation of both graphite and SWNT [|TIJ and correctly 
calculated the elastic modulus, we can compute the strain energy of SWNTs with different 
deformations. Fitting the values of E and Poisson ratio a of SWNT and graphite sheet 
which are calculated from the force-constant model, we obtain the value of So is 57 eV/A . 
Since the bond-length change of SWNT is determined by the first order derivative of energy 
rather than energy itself, although the first term energy of eq. ([7]) in straight SWNT is very 
small m^irjl, the value of Sq still strongly affect the equilibrium structure of SWNT. 



The carbon atom positions of SWNT can be located from the rolling graphene sheet |L4 
the vectors Uij can be easily obtained, 
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r sm6ij)e z - 2p sin —ie rj + posin^e^, (8) 



where is the angle between and the circumference direction. e z , e rj and are the 
unit vectors of axis, radial and tangent direction at atom i, respectively. <pij is the relative 
azimuth between the atom i and the nearest neighbor atom j, independent on the site i, we 
have 

r cos% 

</>H = -■ (9) 

Po 

After considering the relaxation effect, the radius and bond length will be stretched. Uij is 
approximately obtained, if po is replaced by the relaxed radius p, and 4>ij is n °t changed, 
where we have neglected the higher order terms. Then the total strain energy E t can be 
calculated. We can obtain the analytic formula of E t , then minimizing E t calculate e, but the 
obtained formula is very complex. From the numerical calculating, we find the same square 
reverse proportion relation between e and po which is expected by CET. The coefficient is 
slightly dependence on the helicity of tube, since it only depends on an average of cos n 8ij 
in three bonds. Fig. 1 shows the radius dilation of some SWNTs expected by the expanded 
Lenosky model. It is very good in agreement with the characteristics po square reverse 
proportion relation. As comparison, a first principles LDA calculation also be depicted in 
Fig. 1, which is well consistent with our results. Since the structure relaxation can be 
obviously affected by the energy term from bond length change, our obtained value of £, 
from the relaxation is more exact than from the strain energy itself of tube. Therefore, our 
result shows the expanded Lenosky model which £ is fitted from the force constant model 
of graphite can be well used in the studying of the carbon network deformations. 

In our previous paper |TJ| using Tight-Binding (TB) approximation, we have successly 



calculated the strain energy of the bending tube, and simultaneously obtain the Young's 
modulus and the effective wall thickness of SWNT. It clarified the discrepancy on the wall 
thickness of single-layer structure P-LTJ]. The expanded Lenosky model can be used to 
study the same problem, too. Fig. 2 shows the strain energy per atom E b of the (5,5) tube 
as a function of the bending radius R 2 . The data follow quite well the expected behavior 
Eb = E r + A/ R 2 , where E r is the rolling energy of tube. The fitting value of A is about 
185 eVA /atom, slightly larger than the TB result 173 eVA /atom Fig. 3 shows A is 
a linear function of p 2 , a* + b* p 2 , which is consistent with the relation expected by CET, 
\ = Vt Y b (p 2 + & 2 /4), Where Vt = 2.62 A /atom is the occupied area per carbon atom in 
SWNTs, Y is Young's modulus, b is the wall effective thickness, and p is the radius of tube. 
The value of b* is 16.2 eV/atom, slightly larger than the TB results 15.3. Since A is far 
larger than a*, it is difficult to exactly fit the value of a* from the calculated data, due to 
the calculating errors, a* ~ 0.4, if only using the data of (n, n) and (n, 0) tube, but a* « 1.5, 
if using the total data, where the unit of a* is ro 2 eVA /atom . However, our obtained the 
value of b* is insensistive to the selected data, and a* is between 0.4 and 1.8, consist with 



the TB result 1.05 [T^]. Fig. 4 shows the strain energy of straight tube (1/R = 0) E r is 
reverse proportion to the square of tube radius p, the coefficient is about 1.48 eVA /atom, 
good consistent with the TB calculation 1.44 eVA /atom. In previous TB approximation, 
since the unit cell of chiral tube is very large, calculating strain energy is an arduous work, 
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and the method include twelve parameters. We think the expanded Lenosky model can well 
be used to determine the carbon network structure. 

In summary, our studies show the tube relaxation is inverse proportion to the square of 
the tube radius. The expanded Lenosky model which only include four parameters can well 
describe the deformation of the bending carbon nanotube. It may be a simple replacement 
of some complex methods in determining the equilibrium structure and the deformed energy 
of carbon networks. 

The authors acknowledge the discussion with Dr. Haijun Zhou. 
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Fig. 1 

FIG. 1. The radius dilation of SWNT e versus the radius of tubes. The solid line drawn across 
the data corresponds to a least square fit to the behaviors. The a values of the (n, n) (n, 0) and 

(n, 1) tubes are 0.071, 0.059 and 0.062 A , respectively, slightly depend on the helicity of SWNT. 
The results of LAD is from the fig. 1 (c) of ref. Jffl. 
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FIG. 2. Strain energy per atom versus the bending radius R in (5,5) tube. The solid line is a 
fit to the E r + X/R 2 , where E r is the straight (5,5) tube strain energy. A = 92 x 1.42 2 eVA /atom. 
Here tq = 1.42 A, is the bond length of graphite. 
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Fig. 3 

FIG. 3. the value of A of some (n,0) and (n,n) tubes and some chiral tubes. The solid line 
is a fit the all data to a* + b*p 2 , a* « 1.5 r 2 eVA /atom b* ^ 16.2 eV. (TB results: a* ps 1.05, 
b* ~ 15.3 eV in Ref. |jl4|)- The inset shows the fitting solid line from all data nearby p = 0. To 
compare, we give the TB results, too. 
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Fig. 4 



FIG. 4. The strain energy per atom versus the radius of (n, n) and (n, 0) tubes. The solid line 
drawn across the data corresponds to a least square fit to the ^ behavior. C = 0.74 • vq 2 m 1.48 

eVA 2 /atom. C is independent on the helicity of tubes. 
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